arXiv:1502.01468vl [math-ph] 5 Feb 2015 


Brownian Motions with One-Sided Collisions: 
The Stationary Case 

Patrik L. Ferrari* Herbert SpohT Thomas Weiss^ 

February 5th, 2015 


Abstract 

We consider an infinite system of Brownian motions which interact 
through a given Brownian motion being reflected from its left neigh¬ 
bor. Earlier we studied this system for deterministic periodic initial 
configurations. In this contribution we consider initial configurations 
distributed according to a Poisson point process with constant inten¬ 
sity, which makes the process space-time stationary. We prove conver¬ 
gence to the Airy process for stationary the case. As a byproduct we 
obtain a novel representation of the finite-dimensional distributions of 
this process. Our method differs from the one used for the TASEP 
and the KPZ equation by removing the initial step only after the limit 
f —>• oo. This leads to a new universal cross-over process. 
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1 Introduction 


We will study an infinite system of interacting Brownian motions with 
Xn{t) G M, n G Z, denoting the position of the n-th Brownian particle on 
the real line at time t. Initially the positions are ordered as x„(0) < x„+i(0), 
with the convention that a;o(0) < 0 < Xi(0). As indicated in the title, parti¬ 
cle n +1 interacts through a steep, narrowly supported potential with its left 
neighbor, n, only. In the limit of zero support, the singular limit studied in 
our contribution, this interaction amounts to Brownian motion n -|- 1 being 
reflected from Brownian motion n. A mathematical dehnition will be given 
below. Under this dynamics the order is preserved. 


Xn{t) < Xn+l{t) 


( 1 . 1 ) 


for all times t > 0. In previous work we investigate the case of initial condi¬ 
tions with equal spacing, x„(0) = n [H]. Another natural initial condition is 
such to have the process space-time stationary, which is accomplished by as¬ 
suming that {a;n(0), n G Z} is a Poisson point process with uniform intensity 
which, without loss of generality, can be taken as 1. Then {xn{t),n G Z} is 
again an intensity 1 Poisson point process. 

Our interest are the fluctuations of Xn{t) for large t and n. To understand 
their properties one hrst has to hnd out how an initially small perturbation 
close to the origin propagates in time. This path is known as characteristic. 
In our model, because of the one-sided collisions, the characteristic turns out 
to be a straight line with velocity 1. If n = ['dtj, [-J denoting integer part, 
then for'd 7 ^ 1 only the randomness of the initial conditions plays a role and 
the fluctuations of Xn{t) will be Gaussian asymptotically on the scale. 
However close to the characteristic, i.e., n = [t + with r = 0(1), one 

observes non-Gaussian fluctuations in the scale, the properties of which 
will be analysed in great detail in this contribution. With our methods we 
can handle the stochastic process in r at hxed t. Two-time properties along 
the characteristic are known to be difficult. For example, a long-standing 
problem is to obtain the joint distribution of xit](t),xi 2 t]{‘ 2 t) for large t. At 
the time of writing Johansson reports on asymptotic results for the model 
studied in this paper [23] . 

Our results are closely linked to the one-dimensional Kardar-Parisi-Zhang 
(KPZ) equation [21] with stationary initial data. KPZ is a stochastic PDF 
for a height function h{x, t) G M and reads 


d,h = + W 


( 1 , 2 ) 


with W space-time white noise. As written the equation is only formal, but a 
precise mathematical meaning has been given mm- As random initial data 
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h{x, 0) we choose the statistics of two-sided Brownian motion with constant 
drift b. The dynamics is stationary in the sense that x h{x,t) — h{0,t) is 
again two-sided Brownian motion with drift b [16]. Very recently, Borodin 
et al. [8] succeeded in writing down reasonably concise formulas for the dis¬ 
tribution of h{x,t), conhrming the prior replica computation [20]. Through 
an intricate asymptotic analysis they establish (Theorem 2.17 of [H]) that in 
distribution, for hxed r, 

lim (h{—2bt + 2t) + (-jy + b'^)t — 2brt^^^') = ^stat(^"), (1-3) 

where ^stat denotes the Airy process corresponding to stationary initial data. 
In spirit one should think of Xn{t) as h{x, t) with x being a continuum version 
of the discrete particle label n. More precisely, as one of our results we will 
establish in Theorem 12.21 that 

Jhn V^/^(a:Lt+2rt2/3j(t) - 2t - 2rt^/^) = Mstat(’^), (1-4) 

which is the immediate analogue of fll.Op . In fact, convergence is proved 
in the sense of hnite-dimensional distribution, not only for the one-point 
distribution. 

Similar results have been obtained earlier for the stationary PNG model 
[28] and for the stationary TASEP [IT]. For the latter, the full stochastic 
process in r has been worked out [1]. The expression we obtain for the joint 
distributions of Mstat, see Dehnition l2.ll is new and differs from the one in [1]. 

For several reasons we believe that it is of interest to add a third model 
to the list of KPZ type models with stationary initial conditions. Obviously, 
the universality hypothesis is further strengthened. More importantly our 
model provides a bridge to diffusion processes with one-sided interaction as 
discussed in [31]. Besides we also have to develop a method different from 
the previous ones. As in the case of PNG and TASEP, one cannot study the 
stationary initial conditions directly. One has to start from a step, in our 
case meaning that to the right of 0 the Poisson point process has density 
1, while to the left it has density p, p < 1. Surprisingly, as for PNG and 
TASEP, by a Burke type theorem the left-half system can be replaced by a 
boundary condition for xo{t) and, in fact, only the right-half system with 
labels {n > 0} has to be considered. For PNG and TASEP the limit p —)■ 1 
has been accomplished for hxed t, while here we hrst take the limit t —)■ cxo 
at step size 1 — p = This leads us to a novel transition process, see 

Theorem 12.61 The stationary case, 5 = 0, is then reached through a careful 
analytic continuation. 
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2 Main results 


To state our result we introduce the rescaled process 

r ^ Xt{r) = {t) -2t- (2.1) 

and dehne the limit process r Astat(?")- 

Definition 2.1 (Airystat process). Let Pg be the projection operator on [s, cxo) 
and Pg = 1 — Pg the one on {—oo, s). Set 


('® 1 ) ^ 2 ) 


(s2-sir 

e 4(r2-ri) 

^47r(r2 - ri)’ 


and define 


P = 1. — Pg, Vj-i roPso ■ ■ ■ Fr 1 r Ps Pr ni 

*1 *2 > m—lp m 'mpl' 


as well as an operator K with integral kernel 

poo 

K{si, S 2 ) = / dx Ai(r^ -|- Si -1- a;)Ai(r^ + S 2 + x). 

Jo 

Further, define the functions 


( 2 . 2 ) 

(2.3) 

(2.4) 


2^3 


P = Si + eFi / dx dyAfirf + y)e^^'^, 


. 2^3 


./Sl 

bCX) 


f*{s) = —e 3'’i / dx Ai(r^-I-a:)e 

J S 

poo 

g{s) = 1 — / dx Ai(rJ-I-x)e'’^'^. 


(2.5) 


With these definitions, set 

G^{r, s)=n- ((1 - PK)-^ {Pf* + PKPsA + {P- Psfil ), 9 ) , (2.6) 

where {■,■) denotes the inner product on L^(M). Then, the Airystat process, 
Astat, is the process with m-point joint distributions at ri < r 2 < ■ ■ ■ < r^ 
given by 


k=l 


P( p|{Astat(’^fc) < Sk}j = ^ [Gm{r, s) det (1 - PK) 


2 = 1 


^L 2 (M) 


(2.7) 
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We can now state our main result. 

Theorem 2.2. In the sense of finite-dimensional distributions, 

lim Xt{r) = (2.8) 

t^OO 

Remark 2.3. The joint distributions of the Airystat process were first ob¬ 
tained in |1], see Definition 1.1 and Theorem 1.2 therein. In Definition 12. II we 
state an alternative formula for the joint distributions of the Airy^tat process. 
The main difference between the two formulas is that in [1] the joint distri¬ 
butions are given in terms of a Fredholm determinant on T^({1,..., m} x M), 
while here we have a Fredholm determinant on L^(M). A similar twist was 
already visible in [27] and has been generalized in [9]. 

Since {xn(t),n e Z} is a Poisson point process, the process Xfir) — Xt{0) 
is a scaled Poisson jump process up to a linear part and 

lim (Xfir) - XfiO)) = B{2r) . (2.9) 

t^oo ' ' 

Hence the limit process .4,stat(’") — -Astat(O) must also have the statistics of 
two-sided Brownian motion, a property which is not so easily inferred from 
our formulas in Dehnition 12.11 But we will provide a direct proof of this fact 
in Section [HI Note that W(0) and B{2r) are not independent. 

As already familiar from other models in the KPZ universality class imsi 
[T31[Tn] , the proof of Theorem 12.21 proceeds via a sequence of approximating 
initial conditions. Firstly we consider the case where xo(0) = 0 and assume 
that the particles on M+ is a Poisson process with intensity A > 0 and on M_ 
is a Poisson process with intensity p > 0. In other words, a:„(0) = (n, E Z, 
with 

Co = 0, 

Cn - Cn-i ~ exp(A), forn>0, (2.10) 

Cn - Cn-i ~ exp(p), for n < 0. 

As explained in Lemma 17.11 setting (^o = 0 will induce a difference of order 
one as compared to the case considered in Theorem 12.21 In the scaling limits 
such differences are irrelevant. Thus it is enough to prove Theorem 12.21 for 
the initial conditions (12.1 Op with A = 1 = p. In the sequel Xn{t) always 
refers to the initial conditions fl2.10p . in such a way that the choice of the 
parameters A,p can be inferred from the context. We obtain the hxed time 
multi-point distributions of the system {xn(t),n G Nq} in terms of a Fredholm 
determinant in the case A > p. The restriction to non-negative integers comes 
from Burke’s theorem by which the particles with n < 0 can be replaced by 
choosing XQ{t) as a Brownian motion with drift p. 
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Proposition 2.4. Let X> p> 0. For any finite subset S o/Nq, it holds 

pf ^ + det(l-Xa/CXa)L2(5xK), (2.11) 

V„cQ /V P J 


where Xa{^,f) = 1 (^ > o,n)- The kernel K. is given by 

/C(ni,^i;n2,6) = -0 ni,n2 (■Cl) 'C2)l(n2>ni) + Ti{ni, .^i; 77.2, .^ 2 ) 


( 2 . 12 ) 


where 


7,n2(el,6)=p■"^e"^^ 


for 772 > 0, 


(^2 - 

0ni,n2(6,6) = -( -TTrlSl<«2) 1 < JT'l < ^2, 


(2.13) 


(772 - 77i - 1)! 


and 


/C(77i,^1; 772,^2) = 


(27ri)2 


dw ® 


gtto 2 / 2 +^l«; 


dtr 


Pq gt22/2+^22 (^_^)n2 y; _ 2 ;’ 


34 ui 2 / 2 + 5 iui/_„ ,\ni 


iR—£ 


^(^^ 2 , 6 ) = 


27ri 


d^ 


ro,- 


w + \ 

Z + p 


(2.14) 


for any fixed 0 < £ < A. 


Notice that this result holds for X > p only and not for the most inter¬ 
esting case X = p. The latter can be accessed through a careful analytic 
continuation of the formulas. One of the novelty of this paper is to perform 
the analytic continuation after the scaling limit. This allows us to discover a 
new process, called finite-step Airystat process, describing the large time limit 
close to stationarity (actually, one still needs to take care of the random shift 
of Xo(0), which is however irrelevant as it goes to zero after scaling in the large 
time limit). As before, this process is dehned through its hnite-dimensional 
distributions. 


Definition 2.5 (Finite-step Airystat process). The finite-step Airystat process 
with parameter 5 > 0, A[2,t! Is the process with m-point joint distributions at 
ri < r 2 < ■ ■ ■ < Vm given by 


P 




(5) 


atati.rk) < Sfc} = 


k=l 


2 = 1 




(2.15) 
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where Xsi.'f'kix) = l(a; > s^) and the kernel is defined by 


K^{ri, Si;r2, S2) — —K.i,r2(si) 'S 2 )l(ri<r 2 ) + -^ri,r2('Sl; ^2) + ( 5 /ri (Sl)fi'r2 (^2) • 

(2.16) 

Here, K.j,r 2 'Is defined as in /fg.gP . and 


^ri,r2 (®1) ^ 2 ) 


g|r-|+r2S2 roo 


g|r-3+risi 


2^3 


= 1 - e-H^ 


-Trry-risi 


dxe ''^^Ai(r^ + Si + a;)Ai(r2 + S2 + x) 

poo 

/ dx Ai(r^ + Si + 


rkjiu 

(S 2 ) = _ g|r3+r2S2 / Ai(r2 + S 2 + 

Jo 

(2.17) 


As mentioned already above, we are going to take the limit to stationarity 
after the long time limit. However, in general, the limits t —)■ cxo and A —p J, 0 
do not commute. Therefore we have to consider A —p > 0 (to be able to apply 
Proposition 12.4p . but vanishing with a tuned scaling exponent as t —)■ cxo, a 
critical scaling. We set A — p = for 5 > 0. As will be proven with this 

choice the limit t ^ oo commutes with 5 J, 0. 

Such considerations lead naturally to dehne the rescaled process as 


r xl‘\r) = - 2( - 2refi , (2.18) 


where the superscript of x indicates A = 1 and p = 1 — 

The second main result of our paper is the description of the joint distri¬ 
butions of the rescaled process in the long time limit. 

Theorem 2.6. For every 5 > 0, the rescaled process fl2.18p converges to the 
hnite-step Airy^tat process 


lim Wf ^ (r) = Ai2,t {r), (2.19) 

£—)-oo 

in the sense of finite-dimensional distributions. 
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3 Semi-infinite initial conditions 


3.1 Well-definiteness 

Consider the initial conditions stated in fl2.10j) . First we show that the system 
with inhnitely many particles is well-dehned. For that purpose we use the 
Skorokhod representation [21EQ] to dehne the reflected Brownian motions. 
This representation is the following: the process x{t), driven by the Brownian 
motion B{t), starting from a:(0) G M and being reflected at some continuous 
function f{t) with /(O) < x(0) is defined as: 

x{t) = x(0) + B{t) — min {O, inf (a:(0) + B{s) — /(s))| 

= max {x(0) + B{t), sup (/(s) + B(t) — i?(s))}. 

0<s<t 

Let Bn, n E Tj, he independent standard Brownian motions starting at 0 
and define the random variables 

n 

= sup '^{Bi{si+i) - Bi{si)) (3.2) 

with the convention = 0 and = t. We will define the system 

{xn{t),n G Z} as the limit of half-infinite systems {x^T^\t),n > —M} as 
M —>■ oo, where 

= max {Yk,n{t) + Ck], n>-M. (3.3) 

ke[-M,n] 

Notice that these processes indeed satisfy the Skorokhod equation, 

= max {C„ -h Bn{t), sup (xi^i(s) Bn{t) - Bn{s))}, (3.4) 

0<s<t 

for n > —M, while the leftmost process is simply 

= C-M + B_M{t). (3.5) 

Thus as desired xi^'^ (t) is a Brownian motion starting from (n and reflected 
off by x^^l for n > —M. The representation fl3.3p can be sees as a zero- 
temperature case of the O’Connell-Yor semi-discrete directed polymer |2S] 
with appropriate boundary conditions (see discussion at the end of this sec¬ 
tion). 

Next we show strong converge of the system {x^\t),n > —M} to the 
processes we are studying. 






Proposition 3.1. Let us define 


Xn{t) = max{yfc,n(t) + Ck}- (3.6) 

k<n 

Then, for any T > 0, 

lim sup \xl^\t) — Xn{t)\ = 0, a.s. (3.7) 

M^oo jgjO,T] 

as well as 

sup \xn{t)\ < oo, a.s.. (3.8) 

te[o,T] 

For the proof of this proposition we hrst need the following concentration 
inequality, which is Proposition 2.1 of [23] : 

Proposition 3.2. For each T > 0 there exists a constant C > 0 such that 
for all k < m, 5 > 0, 

p/- ^fc,m(T) > 1 + < ^g-(m-A:+l)53/2/c_ 

\2^{m-k + l)T - J- 

Proof of Proposition 17.11 Let us dehne an auxiliary system of processes, 
which we will use later in proving the Burke property, by 

= C-M + B_M{t) + pt, (3.10) 


and 


= max {C„ + Bnit), sup (x[,^J(s) + 5„(l) - B„(s))} (3.11) 

0<s<t 

for n > —M. This system differs from x^\t) just in the drift of the leftmost 
particle, which of course influences all other particles as well (the choice of the 
extra drift is because the system with inhnite many particles in M_ generates 
a drift p). This system of particles satishes 

= max [Y_M,n{t) + C-M, max {Yk,n{t) + Cfc}}, (3.12) 

A;e[-M+l,n] 

with 

n 

YkAt) = sup (psfc+i + ^(Bj(si+i) - Bj(si))y (3.13) 
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Also, we have the inequalities 

hfc,n(t) < Yk,n{t) < Yk,n{t) + pt. (3.14) 

Consider the event 

Am := {Y-M,n{T) + C-M + pT > —M/2} U {Yn^n{T) + Cn < —M/2}. (3.15) 

We can deduce exponential decay of P(Am) in M from combining the expo¬ 
nential tails of Yn^n and (/n with Proposition 13.21 using <5 = 1 and elementary 
inequalities. In particular X]m=i ^’(Am) < cxo, so by Borel-Cantelli, Am 
occurs only hnitely many times almost surely. This means, that a.s. there 
exists a Mt, such that for all M > Mr'- 


Y-M,n{T) -\- Q-m + pT < —M/2 and Yn^niT) -|- Cn > —M/2. (3.16) 

Consequently, Y_M,n iT)+C -M < Yn^n{T) + (n for all M > Mt and therefore 


i„(T) = 

(3.17) 

Moreover. aDDlvinsr f|3.14]l. eives 


Y-M,n{t) + C-M < Y-M,n{t) + C-M < Yn^n{T) + Cn, 

(3.18) 

resulting in 

= x<i'x\T). 

(3.19) 

Repeating the same argument, we see that for every t G [0, T] there exists 
Mt such that Xnif) = = xi!^^\t). Applying Lemma l3.3l then gives 

Xn{t) = xl^^\t) = x^'^\t) for every t G [0,T]. This settles the convergence. 

To see (I3.3[l. which is equivalent to sup^gjo.T] < oo, we apply the 

bound 

n 

\Yk,n{t)\ < V] ( sup Bi{s) - inf Bi{s)) < oo. 

^ ^0<s<t 0<s<t J 

l = k — — 

(3.20) 


□ 

Lemma 3.3. Consider 0 < ti < t 2 and m, Mt^, Mt^ such that 


Xraiti) = x\m'"\ti) = X^m'"\ti), for i = 1,2. 

(3.21) 

ThcTi 

Xm{tl) = xl^"^\ti) = xl^'^\ti). 

(3.22) 


Proof. This is a straightforward generalization of Lemma 3.2 [T5]. □ 


10 






3.2 Burke’s property 

We establish a useful property which will allow us to study our system of 
interacting Brownian motions through a system with a left-most Brownian 
particle. 

Proposition 3.4. For each n < 0, the process 


^n(t') Cn pt 


(3.23) 


is a standard Brownian motion. 

Remark 3.5. Proposition 13.41 allows us to restrict our attention to the 
half-inhnite system. In fact, conditioned on the path of xq, the systems 
of particles {xn(t),n < 0} and {xn(t),n > 0} are independent, as it is 
clear by the dehnition of the system. Then fl3.23l) implies that the law of 
{xn{t),n > 0} is the same as the one obtained replacing the inhnitely many 
particles {xm{t),'rn < 0} with a single Brownian motion Xo(t) which has a 
drift p. This property will be used to derive our starting result. Proposi¬ 
tion 12.41 

Proof of Proposition \‘3.4\ First notice that 

- C-M - pt, (3.24) 

is a Brownian motion. Now assume x^\{t) — (n-i — pt is a Brownian motion. 
By dehnition, 

X^,^\t)-Cn-l = max{Cn-Cn-l+5n(t), SUp {x^^\{s)-(n-lP Bn{t) - Bn{s))} , 

0<s<t 

(3.25) 

which allows us to apply Proposition 13.61 i.e., we have that 

- Cn-l - iCn - Cn-l) - pt = xl^\t) - (n - pt (3.26) 

is a Brownian motion. Since xl!^'^\t) = Xn{t) the proof is completed. □ 

It is clear, that in the case X = p the process fl3.23p is a Brownian 
motion for n > 0, too, i.e., the system is stationary in n. We also have 
stationarity in t, in the sense that for each f > 0 the random variables 
{xnit) — Xn-iit),n G Z} are independent and distributed exponentially with 
parameter p. The following result is a small modihcation of Theorem 2 in j26j . 
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Proposition 3.6 (Burke’s theorem for Brownian motions). Fix p > 0 and 
let Bit), C{t) be standard Brownian motions, as well as ( exp(p), inde¬ 
pendent. Define the process 

D{t) = max + C{t), sup {B{s) + ps + C{t) — C(s))}. (3.27) 

0<s<t 


Then 


D{t) - C- pt 


is distributed as a standard Brownian motion. 


(3.28) 


Proof. Extend the processes B{t), C{t) to two-sided Brownian motions in¬ 
dexed by M. Dehning 

q{t) = sup {B(t) — B(s)C(t) — C(s) — p{t — s)} (3.29) 

— oo<s<t 

and 

d{t) = B{t)q{0) - q{t), (3.30) 

we can apply Theorem 2 [25], i.e., d{t) is a Brownian motion. Now, 

g(0) = sup{—i?(s) — C{s) -\- ps} = sup{\/2i?(s) — ps} = sup {B{s) — -s}, 

s<0 s>0 s>0 2 

(3.31) 

so by Lemma 13.71 it has exponential distribution with parameter p. As it 
is independent of the processes {B{t),C{t),t > 0} we can write g(0) = C- 
Dividing the sup remum into s < 0 and s > 0 we arrive at: 

-d{f) = q{t) - B{t) - q(0) 

= max \c{t) — pt, sup {—B{s) -\- C{t) — C{s) — pit — s)} — cj; 

f 0<s<t 1 

(3.32) 

which is fl3.28p up to a sign flip of B{s). □ 

Lemma 3.7. Fix p > 0 and let B{t) be a standard Brownian motion. Then 

sup(i?(s) — ps) ~ exp(2p). (3.33) 

s>0 

Proof. The random variable 

sup (i?(s) — ps) (3.34) 

0<s<i 

is distributed as a Brownian motion with drift —p started at zero and being 
reflected (upwards) at zero, at time t. As t —)■ oo, this converges to the 
stationary distribution of this process, which is the exponential distribution 
with parameter 2p. □ 
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So Si S2 S 3 t 

Figure 1: A path vr G 11(0, 0;t, 4) (thick black) and the random background 
noise (grey). 


From a stochastic analysis point of view, the system {xn{t),n > 0} satis- 
hes 

Xn{t) = Cn + Bnit) + for U > 1, 

~ (3.35) 

Xo{t) = Bo{t) + pt. 

Here L'^,n> 2, are continuous non-decreasing processes increasing only when 
Xn{t) = Xn-i{t). In fact, L" is twice the semimartingale local time at zero of 
Xn — Xn-i- Notice that i?o(^) is a standard Brownian motion independent of 
{Cn, Bn{t),n > 1}, but not equal to B^it). 

3.3 Last passage percolation interpretation 

One can also view the system {xn(t),n > 0} as a model for last passage 
percolation (or zero-temperature semi-discrete directed polymer). We assign 
random background weights on the set {M+ x Nq} in the following way: 

• White noises dBn on the lines R+ x {n} for n > 1, 

• White noise di?o plus a Lebesgue measure of density p on the line 
R+ X {0}, and 

• Dirac measures of magnitude Cn — Cn-i on (0,n) for n > 1. 

An up-right path in {R_|_ x Nq} is characterized by its jumping points Sj 
and it consists of line segments [s„_i,s„] x {n}, see Figure [H The set of 
up-right paths can then be parameterized by 

n(ti,ni;t2,U2) = {s e < Sn^ < ■ ■ ■ < Sn2 = fo} • (3.36) 
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The percolation time or weight of a path G If is the integral over the 
background weights along the path. Explicitly, we have: 


wiji) — Bq{sq) + pso + ((0 ~ 0-i)lsi-i=o + Bi(si) — (3.37) 


for ni = 0, and for ni > 0, 


“ Ci-i)lsi-i=o + - Bi{si_i)) . (3.38) 


The last passage percolation time is given by the supremum over all paths: 


^{ti,ni)^(t2,n2) ■ sup w(^7l). 


(3.39) 


7fen(ti,ni;t2,n2) 


The supremum is almost surely attained by a unique path if*, called the 
maximizer. It exists because the supremum can be rewritten as a composition 
of a hnite maximum and a supremum of a continuous function over a compact 
set. Uniqueness follows from elementary properties of the Brownian measure. 
Most importantly, from the dehnition, we have 


Xn{t') T(o,0)— >(t,n) • 


(3.40) 


We will use this interpretation in Section [71 however, it also gives some 
connections to different works. Our model can be seen as the semi-continuous 
limit of a more widely studied discrete last passage percolation model (see 
for example [211122]). Moreover, our last passage percolation model is the 
zero temperature limit of a directed polymer model, which has been studied 
thoroughly in the recent past [HEH]. 

For later use we also define a version without boundary weights: 


n 



(3.41) 
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4 


Determinantal structure - Proof of Propo¬ 


sition 12.4 


In order to prove Proposition 12.41 we first start by considering the transition 
probability for a finite number of reflecting Brownian motions with drifts and 
arbitrary initial positions fProposition 14. ip . Then we will set the drift of the 
hrst Brownian motion to p, see Remark l3.51 and we will randomize the initial 
positions fProposition 14.4j) . 


4.1 Transition density for fixed initial positions 

Proposition 14.11 generalizes Proposition 4.1 na, which has been hrst shown 
in [33], to the case of non-zero drifts. 

Proposition 4.1. The transition probability density of N one-sided reflected 
Brownian motions with drift jl from T(0) = C ^ lo x{t) = ^ G at 
time t has a continuous version, which is given as follows. 

P (T{t) E df|x(0) = C) = (4.1) 

where 

n{C,i) = f - CN+i-k,t)], (4.2) 

n=l ' - > - 


and 

-- 

with fi > — min{/ii,. 


27ri Jm+fj, 
■ ,Tn}- 


dw 






(4.3) 


Proof. We follow the proof of Proposition 8 in [33] . The strategy is to show 
that the transition density satishes three equations, the backwards equation, 
boundary condition and initial condition, the latter one being contained in 
Lemma [4.21 These equations are then used for Ito’s formula to prove that it 
indeed is the transition density. 

We start with the backwards equation and boundary condition: 


drt f 1 d‘^ d \ 

7 ^ = 0, whenever Q = Ci-i, 2 < i < N 
dCi 


(4.4) 

(4.5) 
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To see fl4.5p . move the prefactor inside the integral in the {N + 1 — i)-th 
row of the determinant and notice that the differential operator transforms 
Fkj into —Fk+i^i- Consequently, Q = Q-i implies the (A^ + 1 —i)-th being the 
negative of the (7V + 2 — i)-th row. fl4.4p can be obtained by the computation 



Let /: —)■ M be a C°° function, whose support is compact and has a 

distance of at least some e > 0 to the boundary of . Dehne a function 
F : (0, cxo) X —)■ M as 


F{tX)= [ (4.7) 

Jw^ 


The previous identities fl4.5p and fl4.4l) carry over to the function F in the 
form of: 


— = 0, for Ci = 0-1, 2<i <N 
oCi 


(4.8) 

(4.9) 


Our processes satisfy Xn(t) = Cn + f^nt + Bn(t) + where are indepen¬ 

dent Brownian motions, = 0 and L”, n > 2, are continuous non-decreasing 
processes increasing only when Xn{t) = Xn-i{t). In fact, L"’ is twice the semi¬ 
martingale local time at zero of Xn — Xn-i- Now £x some e > 0, T > 0, dehne 
a process F{T + e — t, x{t)) for t G [0, T] and apply Ito’s formula: 


X d 

F {T + e-t,x{t)) = F {T + e,x{d)) + / -—F{T + e-s,x{.s))<ls 


N rt 




a 


. N 


E 


m,n=l ^ 


F {T + e — s, x{s)) da:„(s) 

t q2 

—F (T 4- £ - S, X{s)) d {Xrr,{s),Xn{s)) . 

(4.10) 


From the definition it follows that = /Xndt + d5^(t) + dL’^(t) and 


d ^n(^)) d (^) 7(^) ) ^m,ndt, 


(4.11) 


16 







because continuous functions of finite variation do not contribute to the 
quadratic variation. Inserting the differentials, by fl4.9p the integrals with 
respect to ds integrals cancel, which results in: 


/■* B 

iini) = F (T + £, a;(0)) + V / —F {T + e - s, x{s)) dF^s) 

+ Z / Jipf (r + e-s,f(s))dL”(s). 

n=l'^0 C'U 


(4.12) 


Since the measure dL^(t) is supported on {xn(t) = x„_i(t)}, where the 
spatial derivative of F is zero (see 04.81) h the last term vanishes, too. So 
F (T + e — t,x(t)) is a local martingale and, being bounded, even a true 
martingale. In particular its expectation is constant, i.e.: 

F(r + e, C) = E [F (T + £, x(0))] = E [F (e, x(T))]. (4.13) 

Applying Lemma [4.21 we can take the limit e —?■ 0, leading to 

F(r,C) = B|/(f(T))]. (4.14) 


Because of the assumptions we made on / it is still possible that the distri¬ 
bution of x(T) has positive measure on the boundary. We thus have to show 
that is normalized over the interiour of the Weyl chamber: 

Start by integrating 04. 2 p over G [^Ar_i,oo). Pull the prefactor in¬ 
dexed by n = iV as well as the integration inside the I = 1 column of the 
determinant. The {k,l) entry is then given by: 


g-MivCiv-t/4^/2 / dewe^^«^F„,i(eiv - (N+i-k, t) 


Tiv-i 
,2 




(4.15) 




The contribution from x = ^n-i is a constant multiple of the second column 
and thus cancels out. The remaining terms are zero for k > 2, since all these 
functions F^ 2 have Gaussian decay. The only non-vanishing term comes from 
k = 1 and returns exactly 1 by an elementary residue calculation. 

The determinant can thus be reduced to the index set 2 < k, I < N. 
Successively carrying out the integrations of the remaining variables in the 
same way, we arrive at the claimed normalization. This concludes the proof. 

□ 
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Lemma 4.2. For fixed ( G , the transition density as given by 

04.21) . satisfies 

lim/ n(C,e)/(Ode = /(C) (4.16) 

for any C°° function f: —?■ M, whose support is compact and has a dis¬ 

tance of at least some s > 0 to the boundary ofW^. 

Proof. At first consider the contribution to the determinant in 04.2p coming 
from the diagonal. For k = I the products in 04.3p cancel out, so we are 
left with a simple gaussian density. This contribution is thus given by the 
multidimensional heat kernel, which is well known to converge to the delta 
distribution. The remaining task is to prove that for all other permutations 
the integral vanishes in the limit. 

Let a be such a permutation. Its contribution is 


N 


f if) J_ J_ Fk,a{k) {fiN+l—a{k) Cw+1— 


k=l 


(4.17) 


where we have extended the domain of / to being identically zero outside 
oiWj^. We also omitted the prefactor since it is bounded for ^ in the compact 
domain of /. 

There exist i < j with a{i) <i< crifi). Let 

11^1 = {'f ^ ■Cw+i-crp) - CAf+i-i < -^/2} 

W 2 = {i^ ^v+i-(T(j) — CAf+i-i > ^/2}. 

It is enough to ^tricJAhe area of integration to these two sets, since on the 
complement of Wi U W 2 , we have 


f,N+l-a{i) > Cn+I-i — s/2 > Ca^+I-j “ ^/2 > fN+l-a{j) — (4.19) 


SO we are not inside the support of /. ^ 

We start with the contribution coming from Wi. Notice that by 

FkAU) = (4.20) 

all functions F).; with k > I can be written as iterated derivatives of ^ and 
some exponential functions. For each k ^ i with k > cr(/c) we write F).,(T(fc) 
in this way and then use partial integration to move the exponential factors 
and derivatives onto /. The result is 



d'C f i.f,')Fi^a{i){f,N+l—a(i) CN+l—iyt') .^fc,max{fc,(T(fc)} if,N+l—a[k) CN+l—kyt) 

k^i 

(4.21) 
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for a new C°° function /, which has compact support and is therefore 
bounded, too. We can bound the contribution by first integrating the vari¬ 
ables ^N+i-a{k) with k > a{k), k ^ i, where we have a gaussian factor Fk^. 


I(j4.2ipI ^sup|/(x)| (^w+i— (T(i) Cv+i—nI d^v+i— (T(i) 


'w{ 


n a(k){^N+l-a(k) — CN+l-k,t) \ d^N+l-a{k)- 

k<Ca{k),k^i 

_ (4.22) 

IVI consists of the yet to be integrated ^-components that are contained in 
the set Wi fl supp(/). In particular, W( is compact, so the functions Ffc,cr(fc), 
k ^ i, are bounded uniformly in t by Lemma 14.31 The remaining integral 
gives: 


| fl4.2ip | < const 


r-e/2 




(4.23) 


which converges to 0 as t —)■ 0 by fl4.25p . 

The contribution of 14^2 can be bounded analogously with j playing the 
role of i. The hnal convergence is then given by fl4.24p . □ 


Lemma 4.3. For each e > 0 we have 

poo 

\Fk,i{x,t)\ dx = 0, 
\Fk,i{x,t)\ dx = 0, 


lim 

t^o 


lim 

t^o 


1 <l < k < N, 
l<k,l<N. 


(4.24) 

(4.25) 


In addition, for each 1 < k < I < N the function Fk,i{x,t) is bounded 
uniformly in t on compact sets. 


Proof. Let x < —e, and choose a p which is positive. By a transformation 
of variable we have 


\Fk,i{x,t)\ = 


TTw / dtc e 


/2+xw Wi=l PN+l—i) 

nLl('^ + pN+l-i) 

— f d7; /T^-fe-1 V/2+x4;/W n?=l + VipN+l-i) 


' iR+/i 


(4.26) 


< {2Ti)-^Vt ^ 


|dy|gRehV 2 )^Qy|)^ 




where 5'(|n|) denotes a bound on the fraction part of the integrand, which 
grows at most polynomial in |n|. Convergence of the integral is ensured by 
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the exponential term, so integrating and taking the limit t —)■ 0 gives fl4.25p . 
To see (14.241) . notice that hy I < k the integrand has no poles, so we can shift 
the contour to the right, such that /i is negative, and obtain the convergence 
analogously. 

We are left to prove uniform boundedness of Fk^i on compact sets for 
k < 1. For x < 0 we can use fl4.26p to get 

lm(i,i)l<(2ir)-‘ / (4.27) 

rJ 

for t < 1. In the case x > 0 we shift the contour to negative p, thus obtaining 
contributions from residua as well as from the remaining integral. The latter 
can be bounded as before, while the residua are well-behaved functions, which 
converge uniformly on compact sets. □ 


4.2 Transition density for random initial positions 

To obtain a representation as a signed determinantal point process we have 
to introduce a new measure. This measure P+ coincides with P on the sigma 
algebra which is generated by (k+i—Ck, k E Tj, and the driving Brownian mo¬ 
tions i?fc, fc G Z. But under P_|_, Co is a random variable with an exponential 
distribution instead of being hxed at zero. Formally, P+ = P (8) P^g, with 
P^Q giving rise to (^o ~ exp(A — p), so that P is the result of conditioning P+ 
on the event {Co = 0}. This new measure satishes a determinantal formula 
for the joint distribution at a hxed time. 


Proposition 4.4. Under the modified initial condition specified by P+, 
the joint density of the positions of the reflected Brownian motions 
{xnit), 0 < n < N — 1} is given by 


P+(x(t) e dC) = (A - p)X^ ^e *^V 2 +pCo det [Ffc_i(Cv-z,t)] df (4.28) 


with 


Fk{^,t) 


1 

27ri 



^tw^/2+^w^k 

w + X 


(4.29) 


For the related model, the totally asymmetric simple exclusion process, a 
formula similar to the one of Proposition 14.41 also exists [7] . Here we provide 
a direct proof of it. 


Proof of Proposition \4.4\ The hxed time distribution can be obtained by in¬ 


tegrating the transition density (14. ip over the initial condition. Denote by 
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p+{^) the probability density of x{t), i.e., P+(x(t) G d.^) = p+{^)d^. P+{x) 
equals 


/ dC det [F^Mn-i - (N-k, t)] 

lw^n{Co>o} i<k,i<N^ 


=(A - p)A^-^e-*^'/2+p«o 


dC 


W^n{Co>o} 


X det 
l<k,l<N 




27ri 




=(A - p)A^-^e-*^'/2+^«° 


^ n/. 

ct&Sn k=l 


dta, 


27ri 


- (X^lFp^N-a(k)^kyJ^-^F) 


X dCo... 

Jo JCff-2 


N 


=(A - p)A^-^e-'^'/2+^«° ^ (-l)l‘"l JJ 


dwfc F'^k/‘^e^N-a(k)Wk^^ 


ctGSj^ 


kJi J^R+^k 27ri wi + • • • + Wfc + A 


(4.30) 

Since all Wk are integrated over the same contour, we can replace Wk by tCo-(A:): 
flCTiH = (A - p)A^-^e-*^'/2+p«o 


N 


X 


Ef-i)'"'!! 


(tGSn 


= (A - p)A^-^e-*^'/2+p«o f 

k=i J'^ 


2711 w^(i) H-h Wa(fc) + A 

dtci. 


N 


N 

X E(-i)'”'n 

(tGSn 


_-| J iR+/i 
k 


27ri 


— g*«'l/2g?iV-fc«'fcy^-fc 


(4.31) 


w 


a{k) 


Wad) + • • • + Wa(k) + A 


kSl ^ • • • ^ 

We apply Lemma 14.51 below to the sum and hnally obtain 
p+{x) = (A - TT / ^ det 


k=l 

= (A-p)A'^-'e-‘"'/2+"«° det t) 


Z<Af 


W 


l<k,l<N I 


Wi + X 


(4.32) 

□ 
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Lemma 4.5. Given iV G N, A > 0 and Wi, ..., wn G C \ M_, the following 
identity holds: 


N 


(tSSjv 


W 


-- = det 

U!(7(1) + • • • + UJa-(k) + A l<k,l<N 


k=l ‘"-d) 


w? 


Wi + X 


(4.33) 


Proof. We use induction on N. For = 1 the statement is trivial. For 
arbitrary N, rearrange the left hand side of fl4.33p as 


N 


(tSSjv 


W 


a{k) 


WaiV) + • • • + Wa(^k) + A 


fc=i ‘"-(b 


N 


Y. 

1=1 

N 

Y 


w 


N-l 


—^ E (-1)'^^' TT 

Z=1 + ■ ■ ■ + + A ^^s^^N)=l k=l 

N 


w' 


(k) 


Ula{l) + • • • + UJrj(k) + A 


(fc) 


W 


N 


tci + • • • + wn + A 

«=1 a£SN,a{N)=l 


Y (- 1 )'"' n 


N-l k 

^a{k) 


k=l 


'^a{k) P A 


(4.34) 


where we applied the induction hypothesis to the second sum. Further, 


Wcr(N) + A 


N ^k 


n 


(k) 


ir. + • • • + + A + A 


N 

n 


Wi 


tCi H-h Wat + A tcz + A 


X 


N 


N 


J2 

(tSSjv 


n"’h4+^E(-i)'’'n"’h4 


k=l 


(tGSn 


k=l 


N 

n 


Wl 


Wi-\ -h Wat + A tC; + A 


( det 

k—l+Sk,N 

Wl 

\l<k,l<N 



+ A det [wf ^1 ) , 

T AT ^ '' -1 / 


1<A:,K7V 

with Af being the Kronecker delta. Inserting the identity 
k—l-\-SkN 


det 

l<k,l<N 


we arrive at 


i3®= n 


w, 


N 


= (tCi H-h wm) det \wf ^1 , 

l<k,l<N ‘ 


Wl 


1=1 


, det \wf ^1 = det 
Wl + X J l<fc,Z<Af l<k,l<N 


w' 


wi + X 


(4.36) 


(4.37) 
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To show 04.361) we introduce the variable WAf+i and consider the factor¬ 
ization 


N 


det [wf = TT(M;Ar+i - Wi) det [wf , 
2=1 


(4.38) 


which follows directly from the explicit formula for a Vandermonde determi¬ 
nant. Expanding the determinant on the left hand side along the {N + l)-th 
column gives an explicit expression in terms of monomials in w^+i- Ex¬ 
amining the coefficient of on the left and right hand side respectively 

provides 04.361) . □ 


4.3 Proof of Proposition 12.4 

We can rewrite the measure in Proposition 14.41 in terms of a conditional 
L-ensemble (see Lemma 3.4 of d! reported here as Lemma 14.61) and obtain 
a Fredholm determinant expression for the joint distribution of any subsets 
of particles position. Then it remains to relate the law under P 4 . and P, 
which is the law of the reflected Brownian motions specihed by the initial 
condition fl 2 . 1 Up . This is made using a shift argument, analogue to the one 
used for the polynuclear growth model with external sources um or in the 
totally asymmetric simple exclusion process [41 1141 128]. 


Proof of Proposition 2-4 The proof is divided into two steps. In Step 1 we 
determine the distribution under P_|_ and in Step 2 we extend this result via 
a shift argument to P. 

Step 1. We consider the law of the process under P+ for now. The hrst 
part of the proof is identical to the proof of Proposition 3.5 [15], so it is only 
sketched here. Using repeatedly the identity 


Fk{f,t)=f dxFk+i{x,t), (4.39) 

J —00 


relabeling := fk-ii and introducing new variables ff for 2 < I < k < N, 
we can write 


det 

l<k,l<N 


N+l-l 


F = n (4.40) 

'' ri - - 2<l<k<N 


where V = G M, 2 < / < fc < N\xf < xf_^}. Using the antisymmetry 
of the determinant and encoding the constraint on the integration variables 
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into indicator fnnctions, we obtain that the measure fl4.28p is a marginal of 


N 

const • TT det 


n=2 




det [Fk-i{^j^,t)] 


N 


(4,41) 


const . n , [U-ll?,".*)] 

n=l 


with 


(4.42) 


(j)n{x,y) = l^<y, for n > 2 

Mx,y) = ePy, 

and using the convention that < y always holds. 

The measure 04.411) has the appropiate form for applying Lemma 14.61 
The composition of the (j) functions can be evaluated explicitly as 


i{x,y) = {(l)i*---*(pn){x,y) = 


for n > 1, 


(v — 

0m,n(a;,2/) = (0m+l * ' ' ' * 4>n){x,y) = — -for U > m > 1. 


Define 


K-kiO ■■= 


- 1 ) 


n—k 


{n — m — 1)\ 

ptw'^ 12+^w ^n—k 


27ri 


dtc 


tc + A 


(4.43) 


(4.44) 


for n,k > 1 and some 0 < £ < A. In the case n > k the integrand has 
no poles in the region |t(;| < A, which implies "^^-k ~ (~l)"’~^-^n-fc- The 
straightforward recursion 


(0n*T^,)(o = T:::j_,(o 

eventually leads to condition 04.64p being satished. 

The space 14 is generated by 

{00,n(^?, X),..., 0„_2,n(C-li 0r^-l,n(C"^ ^)}i 

SO a basis for 14 is given by 


(4.45) 


(4.46) 


(4.47) 


Choose functions ^^-k follows 


(- 1 )" 


Z-\-\ 


( C\ ) 27ri ^Fg /2+£,z^n-k + l 

- S (-ir-i r z+A 

I 27ri JFn f,tz^/2+Sz~n-l 


/2+?Z2"-1(2+p) 

24 


2 < k < n, 
k = l. 


(4.48) 



















By residue calculating rules, is a polynomial of order n — k for k >2 

and a linear combination of 1 and for /c = 1, so these functions indeed 
generate 14. To show fl4.66p for £ > 2, we decompose the scalar product as 
follows: 

[ d^n-kiOK-iiO+ [ d^K-kiOK-eiO- (4.49) 

J M— J M.-|_ 

Since n — /c > 0 we are free to choose the sign of e as necessary. For the hrst 
term, we choose e < 0 and the path Fq close enough to zero, such that always 
Re(t(; — z) > 0. Then, we can take the integral over ^ inside and obtain 




(- 1 ) 


k-l 


(2vri) 


dtc ® dz ■ 

To 




(4.50) 

For the second term, we choose e > 0 to obtain Re{w — z) < 0. Then again, 
we can take the integral over ^ inside and arrive at the same expression up 
to a minus sign. The net result of 04.491) is a residue at w = z, which is given 
by 


27ri 


dz z 


= 5, 


kl- 


(4.51) 


ro 


The case £ = 1 uses the same decomposition and requires the choice e > p 
resp. £ < 0, hnally leading to 


fICTD 


27ri 





l-k 


Z + p 


^k,l- 


(4.62) 


Furthermore, both 4>n{^n ^^*4 4>q(.^) are constants, so the kernel has 

a simple form (compare with 04.67^ 1 


n2 

K(n„6;n2,&) = (4.53) 

k=l 

However, the relabeling := included a index shift, so the kernel of 
our system is actually 

/C(ni, 4; ^ 2 ,6) = + 1,4; n 2 + 1,6) 

712 

= -0n4,n,(4,6)l(n,>nO + (^)i(6)• 

k=l 

(4.54) 

Note that we are free to extend the summation over k up to inhnity, since 
the integral expression for $”_^(.^) vanishes for k > n anyway. Taking the 
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sum inside the integrals we can write 

1 


^27ri)2 e*^V2+6^(-^)n2 

(4.55) 


k>l 


3tw)^/2+giui/ 


-W 


.ni 


with 


r]{w,z) = 


z + X ^ z’^-^z + X) 

[w + A)(.s + p) + A) 


(4.56) 


By choosing contours such that \z\ < |t(;|, we can use the formula for a 
geometric series, resulting in 


r]{w,z) = 


z X 


+ 


z X~ X 


{w + X){z + p) {w + X)w 1 — z/w 


+ 


X-p 


(4.67) 


w — z {w + X){z + p) 


Inserting this expression back into fl4.55p gives the kernel fl2.12p . which 
governs the multidimensional distributions of Xn{t) under the measure P+, 
namely 

P+ ( n {Xn{t) < On} ) = det(l - XalCXa)L^{SxR)- (4-58) 

Step 2. The distributions under P and under P+ can be related via the 
following shift argument. Introducing the shorthand 


8{S,d) := p|{x„(t) < a„}. 


(4.59) 


nGS 


we have 


P+(£^(S', a)) = / dxP+(a:o(0) E da;)P+(T(S', a)|xo(0) = x) 

J M-j_ 

= [ dx{X-p)e-^^-P^^F{£{S,d-x)) 


= P(£(^, 5 _ a;)) / dx a - x)) 

0 ./», dx 


= P(£(^,a)) 


- [ W-Tp(£(S,o)) 

K tZs 


= P(£(S,2))- 


1 d 


P4£(S,S)), 




(4.60) 
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Combining the identity 


P(£(S,S))= (l + ^^j^)p4£(S,S)) (4.61) 

with fl4.58p hnishes the proof. □ 

Lemma 4.6 (Corollary of Theorem 4.2 [ID])- Assume we have a signed mea¬ 
sure on {x”, n = 1,..., A^, i = 1,..., n} given in the form, 

1 ^ 

^^n=l 

where are some “virtual” variables and Zjsi is a normalization constant. 
If Zn 7 ^ 0, then the correlation functions are determinantal. 

To write down the kernel we need to introduce some notations. Define 


[0, ni > n2, 

where {a * b){x,y) = J^dza{x,z)b{z,y), and, for 1 < n < N, 

n-jix) := (0^"’^^^ * J = 1,..., at. (4.64) 

Then the functions 

{.#.'"•“>(1;, x),..., x), x)} (4.65) 


are linearly independent and generate the n-dimensional space 14. Define a 
set of functions {^^_j{x),j = 1,... ,n} spanning 14 defined by the orthogo¬ 
nality relations 

[ dx<I>:_,{x)n-j{^) = kj ( 4 . 66 ) 

Jr 

for 1 < i, j < n. 

Further, if x) = c„<ho (x), for some Cn ^ d, n = 1,..., N, then 

the kernel takes the simple form 

n2 

K(n„x,;n„x,) = x,) + J2'‘:U(x,)KU(x2). (4.67) 

k=l 
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5 Asymptotic analysis - Proof of Theorem 12.61 


Remark 5.1. The change in variables 

X X~^X t 1-4 (5.1) 

reproduces the same system with new parameters A = 1 and p = |. We can 
therefore restrict our considerations to A = 1 without loss of generality. 

Fix A = 1 from now on. According to fl2.18p we use the scaled variables 

rii = t + 

e* = 2t + (5.2) 

p=l- 

with (5 > 0. Correspondingly, consider the rescaled (and conjugated) kernel 
/C"®"'=(ri, si; r 2 , Sa) = 77,3, 6 ), ( 5 - 3 ) 


which naturally decomposes into 


/C'-'^^^(ri,Si;r2,S2) = S2)l(n<.d + (^4 "b ^ 2 , Sa). (5.4) 

Remark 5.2. Instead of integrals over Airy functions 02.17p can also be 
written as contour integrals: 




^ri,r2 (®1) ^2} 

/ri(Sl) = 


(27ri)2 

i f 

2tt\ 


dW 


dZ 


/3-\-r2Z‘^ — S2Z 


I g 27ri/3Qgj 


g7ri/3g 


QW3/3+riW^-siW Z — W 


dW- 


e 2-n'i/3oo^ right of 0 


^-{W^/3+riW'^-siW) 

w 


gr2(s2) - 2^. 


-7ri/3g 


dZ- 


'e’^^/^oo, left of S 


^Z^ fS-\-^2Z^ — S2Z 

Z-6 ■ 


(5.5) 

In the integral defining K, the path for W and Z do not have to intersect. 
In addition, the Gaussian part has a representation in terms of an integral 
over Airy functions: 


el ^ 2 +’’2*2 

^l,’'2 (^Ij ^ 2 ) 2 3~i 


dxe '’^^Ai(r^ + Si + x)Ai(r 2 + Sa+ 3^)- (5.6) 


In order to establish the asymptotics of the joint distributions, one needs 
both a pointwise limit of the kernel, as well as uniform bounds to ensure 
convergence of the Fredholm determinant expansion. The hrst time this 
approach was used is in HZ]. These results are contained in the following 
propositions. 













Proposition 5.3. Consider any ri,r 2 in a hounded set and fixed L. Then, 
the kernel converges as 

lim si; rs, Sa) = K\ri, Si; rs, S 2 ) (5.7) 

t—^OO 


uniformly for (si,S 2 ) G [—L,L]^. 

Corollary 5.4. Consider ri,r 2 fixed. For any L there exists to such that for 
t > to the bound 

|/C''®®‘'(ri,Si;r 2 ,S 2 )| < consti, (5.8) 

holds for all (si, S 2 ) G [—L, L]"^■ 

Proposition 5.5. For fixed ri,r 2 ,L and 5 > 0 there exists to > 0 such that 
the estimate 

|/Cr'^(ri,Si;r2,S2)| (5.9) 

holds for any t > to and Si,S 2 > 0. 

Proposition 5.6 (Proposition 5.4 of [I5]). For fixed ri < r 2 there exists 
to > 0 and C > 0 such that 

\fiZrM,S2)\<Ce-\^^-^^\ (5.10) 

Now we can prove the asymptotic theorem: 

Proof of Theorem A2.(k The joint distribntions of the rescaled process x[^\r) 
under the measure P+ are given by the Fredholm determinant with series 
expansion 

^ k=l ' 

(-1)^ m N 

= E / n ’ 

N>0 ■ fc=l “ 

(5.11) 

where n* and fi are given in fl5.2p . Using the change of variables 
CTfc = t“^/^(Cfc — 2t — and a conjugation we obtain 


ira = E^)T- E /ifd'^^* 

Af>0 ■ ii,...,ijv = l k=l 


Ck>Si,) 


X det 
l<k,l<N 


KX^’'{rk,(Tk]ri,ai) 


(1 + ^1) 


2\m+l-ii 


(1 + 


(5.12) 
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where the fraction inside the determinant is the new conjngation, which does 
not change the valne of the determinant. Using Corollary 15.41 and Propo¬ 
sitions 15.51 15.61 we can bonnd the (fc, /)-coefficient inside the determinant 
by 


consti 


(1 + af 

(1 + ’ 


(5.13) 


assnming the are ordered. The bonnds 

(1 + ^ 


[l + y^y 1 + r 


for i < j, 
for j > 1, 


(5.14) 


lead to 

05.131) < const 4 -- K. (5.15) 

1 

Using the Hadamard bonnd on the determinant, the integrand of 05.12p is 
therefore bonnded by 


N 


dak 


constf n 


k=l 


which is integrable. Fnrthermore, 




Ar>0 


m 


(5.16) 


(5.17) 


which is snmmable, since the factorial grows like {N/e)^, i.e., mnch faster 
than the nominator. Dominated convergence thus allows to interchange the 
limit t ^ oo with the integral and the inhnite sum. The pointwise conver¬ 
gence comes from Proposition 15.31 thus 


lim P+ 

t^OQ 


m 

n <sk} 


det (l - XsK\s) 


L2({ri,...,r,„}xK) ‘ 


(5.18) 


It remains to show that the convergence carries over to the measure P. 
The identity 

|i=t-V3^r>(l-p) (5.19) 

leads to 

p( n s I’d) = (i+J E p* (n s 3d). 

1 2—1 k —1 

(5.20) 
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Notice that in (15.12p . Sj appears only in the indicator fnnction, so differenti¬ 
ation just results in one of the ak not being integrated but instead being set 
to Si- Using the same bounds as before we can again show interchangeability 
of the limit t —?■ oo with the remaining integrals and the inhnite sum. □ 

Before showing Propositions 15.31 and 15.51 we introduce some auxiliary 
functions and establish asymptotic results for them. 


Definition 5.7. Define 


at{r,s) ■= — [ 

2vri JiK 

27ri JiK 

271-1 iro 

= / ^^g-i(z2_l)/2-(2t+2l2/3r+il/3s)(^+l)^_^^-i-2t2/3r 

27ri Ito 

Lemma 5.8. For fixed r and L the limits 

a{r, s) := lim at{r, s) = Ai(r^ -f 

t^OO 

(3{r, s) := lim f3t{r, s) = —Ai(r^ + 


(5.21) 


(5.22) 


t^OO 


hold uniformly for s G [—L, L\. 

Proof. Let H„(a;) be the normalized n-th order Hermite polynomial with 
respect to the weight e“* i.e., satisfying 


H^(x)H„(x)e-"'/2da; = S„ 


(5.23) 


There exist two explicit integral representations for these polynomials: 

1 


Hn(x) = 

HAx) = 


i(27r)T4y^ 

VnJ 

i(27r)®U 


/ dw 

J iR+e: 


(5.24) 


ro 


From [52], P- 201, one obtains the asymptotic behaviour 

^i/i2g-x2/4H^ (x) = Ai(M) -I- X = 2y/n + (5.25) 
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uniformly for any bounded u. 

Introducing the shorthands x = n = t + and 

applying the change of variables w we can write 

7 - 1/3 

(5.26) 

Using Stirling’s approximation and Taylor expansion in the exponents one 
can further analyze this as 

a,(r, s) = ^-lr^-rs+o{t-y^) ^ c>(t-i/3)) (5.27) 

with with the error terms being uniform for s G [—T, L\. The observation 

u = 'n}^^{x — 2^/n) —)> + s, (5.28) 


as t ^ oo, settles the convergence of at- 

Using the second integral representation of the Hermite polynomials one 
can rewrite (3t, too: 

^,(r, s) = -t^/3^27r)^/^e*/2-i'^"^r/2(n!)-i/2H„(a;). (5.29) 

Analyzing the prefactor as done before hnishes the proof. □ 


Lemma 5.9. For fixed r, there exist to, L such that for all t > to and s > L 
the following bounds hold 


|Q:t(r,s)| < e ^ 
\fit{r,s)\ < e"" 


(5.30) 


Proof. We start by analyzing fit- Dehning functions as 


- l)/2 - 2(z + 1) - In(-z) 

f 2 {z) =-2r{z + 1+ \n{-z)) (5.31) 

fi{z) = -s{z + 1), 


we can write G{z) = tf 3 {z) + t‘^^^f 2 {z) + fi^^ffiz), leading to 


A(r, s) 


t^ 

27ri 



dze^^^\ 


(5.32) 


Dehne a new parameter oo given by 

cj = min e} , (5.33) 
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Figure 2: The contour F = 71 U 72 (-R) U 72 (-R) U 73 (i?) used for obtaining the 
uniform bounds. 


for some small, positive e chosen in the following, and let 6 G ( 7 r/ 6 , 7 r/ 4 ). As 
shown in Figure |21 we change the contour Fq to 71 U 'J 2 {R) U 'J 2 {R) U 
with 

7 i = { —l + a;(l + iutan^),M G [—1,1]} 

72 (i?) = { — 1 + u G [oj/ cos^, R]} (5.34) 

^s{R) = {-i + Re^\ue [-e,e]}. 

Since we will only estimate the absolute value of the integrals, the direction 
of integration does not matter. 

If t and s are hxed, the integrand is dominated by the exp (— 2 ;^) term 
for large \z\. Thus the contribution coming from 73 (i?) converges to 0 as 
R —)■ 00 . With 72 = lim/j^oo 72 (-R) our choice for the contour of integration 
is now 7 i U 72 U 

We start by analyzing 


tV3 

27ri 



= eG(zo) 




[—u) tan tan 0] 


du QG{zo+iu)-G{zo) 


(5.35) 


where Zq = —1 + oj. 

Let us consider the prefactor at hrst. Since oj is small we can use 

Taylor expansion, as well as fl5.33lh to obtain the bounds 


f/ 3 (^o) = t (^u^/3 + 0{u:'^)) < (1 + 0{e)) 

fRf.izo) = (a;2 + 0[u?)) < (1 + 0[e)) 


All error terms are to be understood uniformly in s,t,r. The fi term dom¬ 
inates both / 2 , if L is chosen large enough, and f^, for e being small. This 
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results in 


gGizo) ^ ^ ^ 


(5.37) 


To show convergence of the integral part of fl5.35p we hrst bound the real 
part of the exponent: 

Re [G[z(j + in) — G[zq)) 


=Re 

=t 


V? — 2zom , Zn + 9/q ^ f. . Zo + iw\ 1/Q . 

-^- 2 m - In . 2 r in + In - - t^'^siu 


zq 


Zq 


u 


- - In 1 + 


<4 ( 1 - 4 + 


2 2 
2 / 1 


U 


+ In ( 1 + 


u 


u 


Zq 2zq 


u 


u 


7] satishes: 


ti/3 


r] = 


u 


{i-uy 


2(l-n;)V (l-a;)2 


= t'/3^(l + C>(a;))-r(l + C>(a;)), 


(5.38) 


(5.39) 


where we used |n| < u. Given any e we can now choose both L and to large, 
such that the hrst term dominates. Consequently r] will be bounded from 
below by some positive constant rjo. The integral contribution coming from 
7 i can thus be bounded as 


IdSSSDI =e 


_ ^GPo) 


< e 






[—u tan 9^u) tan 6] 


r ^ 1 , 3 / 2 ^ r ( 5 . 40 ) 


271 


271 


-1,3/2 1 

= e 


< e 


2v^ 

Finally we need a corresponding bound on the 72 contribution to the 
integral. By symmetry this case covers also the contour 7 ^. Write 


27ri 




' 72 


27ri 


id 


du eGbiW«)-G(.i)^ 


(5.41) 


with zi = —1 + 0/(1 + itan 6 *). From the previous estimates one easily gets 


|gG(zi)| < gG(zo) < 


(5.42) 
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so the remaining task is to show boundedness of the integral part of fl5.4ip . 

At hrst notice that the real part of the /i contribution in the exponent is 
negative, so we can omit it, avoiding the problem of large s. By elementary 
calculus, we have for all u > a;/ cos6', 

^Re(/3(-l+W^)) <0, (5.43) 

that is, 72 is a steep descent curve for /s. We can therefore restrict the contour 
to a neighbourhood of the critical point zi, which we choose of magnitude 5, 
at the expense of an error of order 




du gG'bi+“'=‘®)-G'bi) 


< (1 + / du 


7/3(Me‘®)+t2/3/2(ue'®) 


(5.44) 


where fi{z) = fi{zi + z) — fi{zi). Taylor expanding these functions leads to 


Re(f/ 3 (ue‘')) = fRe(e3'^)u^ (1 + 0{S)) (1 + 0{e)) 

cos"' o 

Re(f2/3/2(Me‘^)) = 2Re(e2‘^)f2/3ru^ (1 + 0(6)) (1 + Oie)) 

cos 6 

< X2k| ■ 


for some positive constants X 2 : X 3 i by choosing 6 and e small enough. For 
large L and to, dominates over X 2 \r\, so we can further estimate: 


[ du 

gt/3(ue'®)+l2/3/2(ue‘®)+/o(ue*®) 

Jo 






2 

Xstu'^' 


(5.46) 


This settles the uniform exponential bound on /3t. 
Combining fl5.26p and fl5.29p one obtains 


Q'f (^, _ r (^\„—2rs—s^t~^^^ 

n ( \ Jty 

Pt{r,s) 


(5.47) 


for some function ft{r). From the convergence of at and (3t h is clear that 
ft converges, too. Since we already know that (5t is uniformly bounded by a 
constant times the exponential bound on at follows. □ 


Proof of Proposition 15. ^ Regarding the hrst part of the kernel, we notice, 
that Uj = 0 does not appear in our scaling, so we can use the formula (for 

77,2 > ni) 


4^ni,n2 ('bl, '^ 2 ) (0ni * . . . 0ri2—l) ('Cl, 'C 2 ) 


(C 2 - 

(u2-r7-l)! 


(5.48) 
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This is the same function as in [15], proof of Proposition 5.1, so the limit 


lim 

t^OO 


ri,r2 


(si, S 2 ) — 


^47r(^ 


rij 




is not proven here. 

The different parts of the remaining kernel can be rewritten as integrals 
over the previously dehned functions a and /3. For /C, choose the contours in 
such a way that Re( 2 ; — tc) >0 is ensured. 


^i/3g6 


(2vri)2 , 

■00 


(2vri; 


f f /2+ii{w+l) 

Im-e fro e*^V2+6b+i) (-^)^2 yj - ^ 

f f gi(«'^-l)/2+^l(«'+l) ryj\ni roo 

Im-e fro e*P'-i)/ 2 +? 2 ( 2 +i) (_^) 


712 




= - / dxat{ri,si + x)f3t{r2,S2 + x) 

Jo 

Also rewrite f as follows: 


(5.50) 


e = 


27ri 


dta 


J{w^-l)/2+^l(w+l) I 


-W 


\ni 


1 

27ri J-,^ 
1 


i—e 

.,2 


ta + 1 


= 1 + / dta 

«-£-l 


gt(«;2-l)/2+^l(«;+!) 

ta + 1 


(5.51) 


27ri .A, ^ ^ ' 


^ j 

Similarly, 


I—e —1 
dx Q;t(ri, Si + x). 


^"^(n2,6) = Res_9_p + / dx/lt(r2, S 2 + x)e^"^. 


with 

Res_,,_p = _^-l/ 35 )-n 2 _ 

The residuum satishes the limit 


lim Res„ 

t^oo ^ 


_ g5®/3+r2(S^-S25 


(5.52) 

(5.53) 

(5.54) 
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uniformly in S 2 - The prefactor of the last part of the kernel is simply 
^1/3(1 _ p) = Combining all these equations gives 


/Co®®''(ri, si;r 2 , S 2 ) = - / dx at(ri, si + x)/3t(r2, S 2 + x) 


^Sx 


+ (5 1^1 - y dxat(ri,si +x)j ^Res^ _p + J dx/3t{r2, S 2 + x)e 

(5.55) 

Using the previous lemmas we can deduce compact convergence of the kernel. 
Indeed (omitting the r-dependence for greater clarity) we can write: 


sup 

si,S 2 e[-L,L] 


dx at{si + x)l3t{s2 + x) — / dx a{si + x)l3{s2 + x) 


< dx sup \at{si + x)Ptis 2 + x) - a{si + x)(3{s 2 + x)\ . 

Jo 

(5.56) 

By Lemma 15.81 the integrand converges to zero for every x > 0. Using 
Lemma ES] we can bound it by const • e~^^, thus ensuring that 05.561) goes to 
zero, i.e., /C converges compactly. In the same way we can show the conver¬ 
gence of f and Applying the limit in 05.551) and inserting the expressions 
for a and (3 hnishes the proof. □ 

Proof of Proposition 15.51 Since the convergence 05.54p is uniform in S 2 we 
can deduce 

|Res_^_p| < consti • (5.57) 

Inserting this as well as the bounds from Lemma 15.91 into 05.55p results in 

/ CX) / POO \ 

+ 5 ( 1 + j dxe-(*i+"M 


^consti • e ^2*^-f J dxe 

g-(si+s2) _|_ (J ^consti ■ ^ 


< const • 


(5.58) 

□ 
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6 Path-integral style formula 

Using the results from [9] we can transform the formula for the mul¬ 
tidimensional probability distribution of the hnite-step Airystat process 
from the current form involving a Fredholm determinant over the space 
L^({ri,..., Tm} X M) into a path-integral style form, where the Fredholm 
determinant is over the simpler space The result of [9] can not be 

applied at the stage of hnite time as one of the assumption is not satished. 

Proposition 6.1. For any parameters Xk 1 < k < m, satisfying 


0 < Xm < ■ ■ ■ < X2 < Xi < max {rj - ri, 5} , 

Kj 

define the multiplication operator {MrJ){x) = mrfix)f{x), with 


( 6 . 1 ) 


mr,(x) = 


e for X > 0 
for X < 0. 


( 6 , 2 ) 


Writing K^,{x,y) := K^{ri,x;ri,y), the finite-dimensional distributions of 
the finite-step Azr^stat process are given by 

/ m 

p( n < Sk} 


k=l 


1 d 

^ i=l * 


1,2 (R) 


with 


Q — -\- Pg, Vj., ro-Pso ' ' 'Vr T r Ps n Pf i 

ri ' *1 'lp2 *2 'm—li'm 'mpl Tp * 

where P* = 1 — Pj denotes the projection operator on (—oo, s). 


(6.3) 

(6.4) 


Remark 6.2. The operator Vr^^n for U < Vj is dehned only on the range of 
K^. and acts on it in the following way: 


^rj,riKri,rk ~ Prj,rki ^rj,rifri — frj- 


(6.5) 


In particular, we have also 14- ,^1 = 1- 


Proof. We will denote conjugations by the operator M by a hat in the fol¬ 
lowing way: 

fn = MrJn, 

9n = (6.6) 


Vr-r- = 




Pri,rj = Mr^Kr^,rjM^. 


-1 
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Applying the conjugation also in the determinant in fl2.15p . the identity we 
have to show is: 


det 1 - XsK Xs 




= det (1 - Kt + Ps^Vr^,r 2 Ps 2 • • • K 


ri 


P V 


(6,7) 


L2(R) 


This is done by applying Theorem 1.1 [5]. 

It has three groups of assumptions we have to prove. We merged them 
into two by choosing the multiplication operators of Assumption 3 to be the 
identity. 


Assumption 1 

(i) The operators Ps^Kn, and for < Vj 

preserve L^(R) and are trace class in L^( 


(ii) The operator Wi.nPri 

serves L^(M) and is trace class in L^(M). 


P„.I4-. r ^.Pq ■ K- 1 r Pq K- prC- 

' zi'z-\-l *1+1 'm — ly'm *m 'mpl Tp ir 


Assumption 2 

(i) Right-invertibility: 

(ii) Semigroup property: = K-i,rfc 

(iii) Reversibility relation: Vr^^rjPrj = PrPn^rj 

The semigroup property is clear. To see the reversibility relation, start 
from the contour integral representation fl5.5p of ,r and and use the 
Gaussian identity: 






^-nW'^+xW 


( 6 . 8 ) 


This results in Vj. r Kt 

ii j I j 


Kri,rj + G) On the other hand we have 


dz 


^y4:7^{rj - Ti) 


Q-i^-y)‘^mrj-ri)^riZ'^-zZ _ ^rjZ'^-yZ 


(6.9) 


SO KtVn^rj = Kj-.^rj + which proves Assumption 2 (iii). Noticing 

Remark 16.21 the right-invertibility follows immediately. 
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Assumption 1 (ii) can be deduced from Assumption 1 (i) as shown in 
Remark 3.2, [5]. Using the previous identities we thus are left to show that 
the three operators for A < rj, as well as Ps^Kr^^r and Psifn 

for arbitrary r*, Tj G M, are all L^-bounded and trace class. 

First notice that Vr^^rj{.x,y) = Vo^rj-ri{.—x,—y). Using the shorthand 
r = Tj — Ti and inserting this into the integral representation fl5.6p of V we 
have 


Vri,rj{x,y) = f d\ Ai{—x + X)e^^ 

Jr 


with the new operators 


y + A) = (V'V®) (x.y), 
( 6 . 10 ) 


{x, y) = Ai(—X + y) 

y) = Ai{r‘^ + x — y). 


Introducing yet another operator, {Nf){x) = exp (—(xi + Xj)a;/2)/(x), we 
can write 


( 6 . 12 ) 


The Hilbert-Schmidt norm of the hrst factor is given by 

[ dxdy\{PsMry^^'>N-^){x,y)\' 

Jr^ 

poo p 

— / dx / dymUx)Ai^(—X + 


(6.13) 


'Si 


'Si 


dxm^.(x)e*^^‘'''^^'^'^ / dz . 


The asymptotic behaviour of the Airy function and the inequalities 
Xi > Xj > 0 imply that both integrals are hnite. Similarly, 


dx dy 


(iVU;2)M-i)(x,|/) 

= 6 ^'’^ [ dxdye~^^^~^^^^^e'^^^^~'^'^Ai^{r'^ + x — y)m~‘^{y) 


= 63’’ / dze + z) / dym^^{y)e Ai+Xj)y ^ 

Jr Jr 

(6.14) 

where we used 2r > Xi + Xj cls well. As a product of two Hilbert-Schmidt 
operators, Ps^Vr^^rj is thus L^-bounded and trace class. 
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(6.15) 


We decompose the operator Kr^^r as 

P..K„,r, = (P..Mr,K'_„PMPoK'M-'] 


where 


K^x, y) = Ai{r‘^ + x + y). 

Again, we bound the Hilbert-Schmidt norms, 

[ dxdy \{Ps^Mr^K'_^.Po){x,y)\‘^ 


(6.16) 


_ 4„3 

= e 30 


/ --r? 

< e 30 


dx / dyml.{x)e Ai^{rf + x + y) (6.17) 


'Si 


dxm^.{x) I dze '^'^^^Ai^{r‘f + z)<oo, 


as well as 

/ da: dy 


{PoKM-^){x,y) 

= f da: f dy + x + y)m~^{y) (6.18) 


= / dym^‘^{y) I dz + z). 

Jr ' Jy 

The superexponential decay of the Airy function implies that for every 
Cl > \rj\ we can hnd C 2 such that e^''^'^Ai^(r| + z) < C2e~^^^ . This proves 
hniteness of the integrals. 

Regarding the last operator, start by decomposing it as 

PsJn ® drj = {PsJn «) 0)(0 «) drj) (6.19) 

for some function 0 with L^-norm 1. Next, notice that 


dx dy I {PsMrJn ® (t)){x,y)\^ = / da:(a:)/^ (a:) 


( 6 . 20 ) 


• Si 


It is easy to see that Umg^oo fnis) = 1, so fn is bounded on the area of 
integration. But then the term ensures the decay, implying that the 
integral is hnite. Furthermore, 


dx dy 


(<I>(S)(x,y) = / dym^. (y)g^.(y). 


( 6 . 21 ) 


Analyzing the asymptotic behaviour of grj we see that for large positive 
arguments, the hrst part decays exponentially with rate —6 and the second 
part even superexponentially. 6 > Xj thus gives convergence on the positive 
half-line. For negative arguments, it is sufficient to see that Qrj does not grow 
faster than exponentially. □ 
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7 Analytic continuation - Proof of Theo¬ 
rem 12.21 


First of all let us show that the choice of Xo(0) = 0 is asymptotically irrele¬ 
vant. Denote by the rescaled process as in fl2.18p . where xo(0) = 0, 

and Xt{r) the rescaled process as in fl2.ll) . where —a;o(0) ~ exp(l). This 
corresponds to a hnite shift of the system, which is therefore irrelevant in the 
large time limit. 

Lemma 7.1. For any m G N, ri < r 2 < ... < and Si,..., Sm, it holds 


hm 1P( {Xt{rk) < Sfc} j = hm P n < .4 ) (7.1) 

''■fc=l ' ^k=l 


Proof. We can construct the processes Xn'^ and Xn on the same probability 
space so that, for any n G Z, Xn\t) = Xn{t) — xo(0) and with a:o(0) being 
independent of Xn{t) — xo(0). After scaling we have xj^\r) = Xt(r) — 
xo(0)t~^t^. As converges to 0 is distribution, the result follows. □ 

We know from Theorem 12.61 and Proposition 16.II that: 

( m \ / 1 H \ 

n {^'’(>-4 <.4 = 1 + 53 det(l - PK‘J. (7.2) 

k=l ' ^ i=l 

In this section we prove the main Theorem 12.21 by extending this equation 
to (5 = 0. The right hand side can actually be analytically continued for all 
5 G M (see Proposition 17.4p . Additionally we have to show that the left hand 
side is continuous at 5 = 0. This proof relies mainly on Proposition 17.21 
which gives a bound on the exit point of the maximizing path from the lower 
boundary in the last passage percolation model. 

Proof of Theorem l^.gt We adopt the point of view of last passage percolation 
discussed in Section 13.31 The superscripts of x, L and w indicate the choice 
of p, while A is always hxed at 1. It is clear that for any path vf the weight 
is non-decreasing in p. But then the supremum is non-decreasing, 
too, and: 

x^^\t) <x\^\t), (7.3) 

for p < 1. We know that there exists a unique maximizing path 
if* G 11(0,0;f;77,). We can therefore dehne Znit) := Sq, the exit point from 
the lower boundary specihcally with p = 1. We want to derive the inequality 

< xif\t) + (1 - p)Zn{t). (7.4) 
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(7.5) 


This can be seen as follows: 

^(oloH(t,n) - (1 - p)Zn{t) = ^ snp - (1 - p)Zn{t) 

7rG^(0,0;^,n) 

= ^*-^^( 7 ?*) — (1 — p)sq = ^*-^^( 7 ?*). 

Note that tt* maximizes tc*'^^(7f) and not necessarily In particular 

we have 

< sup (7.6) 

5fen(0,0;7,n) 

Combining the last two equations results in 07.41) . 

07.31) and 07.4p imply that for the rescaled processes xj^\ see 02.181) . we 
have 

X<‘>(r) < xP(r) < xP(r) + 

For any e > 0 it holds 


(7,7) 


( m \ ^ rn \ 

< sfc}) > p( ^ 

k=l ^ k=l ' 

> p( niAf (rfc) + 5t 2/^Zi+272/3r(^) < Sfc}) 

^k=l ^ 

/ m \ m 

> p( ntN'Fi) (*-"/’Z,+j,v.,(«) > e) 


(7.8) 


k=l 


k=l 


Then, taking t —)■ cx), we obtain 


p( nizis. (rk) < Sfc}^ > limsupP (n{A,'“' (rk) < Sfc}^ 

kfc=i / Vfc=i / 

( m ^ 

P|{xfVfc) < Sfc} 

k=l ' 

( m \ in 

(rk) < Sk- e} \ - ^limsupP (Z^^2t'^/3^(t) > t'^^^e/S) . 

k=i P fc=i 

(7.9) 

Using Proposition 17.21 on the last term and Proposition 17.41 on the other 
terms, we can now take the limit <5 —?■ 0, resulting in 

p(nM stat (rk) < Sfc}^ > limsupP (n{.Y<“> (rk) < Sfc}^ 

> l.mmfp( niAf (rk) < > P (flM stat (rk) <Sk- £}"). 

°° k=l ^ k=l ^ 


(7.10) 
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Continuity of fl2.7p in the Sk finishes the proof. □ 

Proposition 7.2. For any r G M, 

lim limsupP (.^t+2t2/3r(^) > = 0. (7-11) 

Proof. By scaling of t and f3, 07.111) is equivalent to 

lim lim sup P (Zt{t + 2t^l^r) > = 0, (7.12) 

t—>oo 

for any r G M, which is the limit we are showing. We introduce some new 
events: 

My := {Zt{t + 2t^l^r) > 

Ey := {L(o,o)^(/3t2/3,o) + -^(^f2/3,o)^(t+2t2/3r,p < 2t + 2f^/^r + (7.13) 

Ny := {.h(o,o)^(t+2t2/3r,t) < 2t + 2t^/^r + 

Notice that if My occurs, then 

-^(0,0)^(t+2t2/3r,i) = -^{0,0)^-(/3t2/3,0) + -^(/3t2/3^o)^-(t+2t2/3r,t)) (7-14) 

resulting in My F Ey <Z Ny. We arrive at the inequality: 

P(M^) = T{My n Ey) + P(M^ n E'pj < P(iV^) + P(E^). (7.15) 

We further dehne new random variables 
At) ^ 

‘^spiked 

fW _ 

SGUE “ 

_ 

SN “ 

By Theorem 7 [3l], for any hxed r G M, 

^gue-^^gue, (7.17) 

where ^gue has the CUE Tracy-Widom distribution, ^^pj^ed follows the distri¬ 
bution of the largest eigenvalue of a critically spiked CUE matrix, as will be 
shown in Lemma 17.31 has the distribution of a standard normal random 
variable .^n for any /? > 0, t > 0. 


L(^42/3^0)^-p+2t2/3,.,q -2t- 2t^/^{r - (3/2) 


T step 

'^(0,l)^>(t+2t2/3r,p 


tV3 

- 2t - 2f2/3r- 


+ (r-/3/2)^ 


fl/3 


+ r^. 


(7.16) 


qo,o)^(/3t2/3,o) 


- /3f2/3 
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Combining these definitions, we have: 


= lP(v^«g’ + SSw S (»■ - fiP? + «)■ 


(7.18) 


Fix s = 3r^ — /S^/16, such that: 




r -+ s = 4r — rB -\ -1-> 

2 ^ 16 8 - 


(7.19) 


Using the independence of and obtain 

¥{Ef,) > p + eSked ^ + f^) ^ p ^ ^ 


= p ( i'S < 


^3/2 

16 


p I Ciw S ^ 


Further, the inequality 


rStep ^ j 

'^(0,l)^p+2t2/3r,t) — -^(0,0)^(t+2t2/3r,p 


leads to 


P(A'^)<P(«gh<4r'-/3Vl6). 

Inserting fl7.20p and fl7.22p into fl7.15p . we arrive at 

^3/2 


P(Ms) < P < 4r= - Q + 1 - P (sg’ < 
By 07.17p and Lemma 17.31 we can take limits: 


16 


^spiked — 2g 


(7.20) 

(7.21) 

(7.22) 

0' 


^ ( dpiked — ]^g 


(7.23) 


0 < limsup limsup P (Mg) 

13^00 t^oo 


< lim 


P ^GUE < 4r" 


0 

16 


+ 1 - P ( < ^ ) P ( espiked(/9) < 


16 


0 

16 


= 0 . 


(7.24) 

□ 


Lemma 7.3. Let r G M 6e fixed. For any fB > 2(r + 1), as t —)■ oo, the 
random variable 

eL. = + (. - /l/2)^ (7,25) 


tl/3 


45 


















converges in distribution, 


Cw ^ f.pik=d(/3). (7.26) 

In addition, ^spiked(/9) satisfies 

lim P (espiked(/5) < fiVlQ) = 1. (7.27) 

p^oo 

Proof. The family of processes T(^j2/3 o)^(/3t2/3_,_4 indexed by n G No and 
time parameter t > 0 is precisely a marginal of Warren’s process with drifts, 
starting at zero, as dehned in [13]. In our case only the hrst particle has a drift 
of 1, and all the others zero. By Theorem 2 [T3], the hxed time distribution of 
this process is given by the distribution of the largest eigenvalue of a spiked 
n X n GUE matrix, where the spikes are given by the drifts. 

Thus we can apply the results on spiked random matrices, more concretely 
we want to apply Theorem 1.1 |6|, with the potential V{x) = —x‘^/2. Since 


L* ■— 7k(;3i2/3_0)^-(t+2t2/3r,n) (7.28) 

represents a n x n GUE matrix diffusion M(t) at time t = tp 2t^/'^{r — 13/2), 
it is distributed according to the density 


Pn{M) = — exp 


Tr(M-flii)^ 

2f 


(7.29) 


where In is a n x n matrix with a one at entry (1,1) and zeros elsewhere. 
In order to apply the theorem we need the density given in equation (1) [6], 
i.e., consider the scaled quantity L*/\/ni. The size of the hrst-order spike is 
then: 

a = t/V^ = sjl + 2t-y^{r - 13/2) = 1 + (r - [3/2)t-^/^ + ^(r^/^). ( 7 . 30 ) 
We are thus in the neighbourhood of the critical value ac = 1. For a > 0, let 

CPf)= f e^^Ai{x + f)dx. (7.31) 

J —00 


With To(s) being the cumulative distribution function of the GUE Tracy- 
Widom distribution, and J’7o,o(si, S 2 ) as in (12.17p . dehne: 

Ffis- a) = Fo{s) (1 - ((1 - PsKo,oPs)-^Ca, P.Ai) ). (7.32) 

Applying (28) [6], we have 

n^l\L*/^t - 2) ^ espiked(/5), (7.33) 
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with 


(7.34) 


P (6piked(/5) < /9Vl6) = Fi(/3Vl6,a), 
where a = /9/2 — r. Since in onr case a > 1, we can estimate: 

|C„(OI < (7.35) 

J-co « - 1 

Combining this with the usnal bonnds on the Airy kernel and the Airy fnnc- 
tion, we see that as /9 —)■ oo, the scalar prodnct in fl7.32p converges to zero 
and we are left with the limit of Fq which is one. 

On the other hand, 

/ y/ni — 2) < s L* < \/?7i(2 + ?7,“^/^s), (7.36) 


and 

V^{2 + n-^/^s) = 2t + -13/2) +(s - (r - 13/2f) + 0(1), (7.37) 

from which the claim follows. □ 

Proposition 7.4. The function 5 ^ (5“^det(l — VK/f) can he extended 
analytically in the domain <5 G M. Its value at 5 = t) is given by 

Gmir, s) det (1 - VK) ^ 2 (r) . (7.38) 

Proof. We nse the identity det(l + A) det(l + B) = det(l + A + B + AB) 
and Lemma [7.51 to factorize 

(5-^ det (1 - VK/J = 6-^ det(l - VK/J = det (l - VK - 6Vfr, <8) 

= det (1 - 5(1 - VK)-^Vfr, ® Qr,) ■ det (l - VK) 
= (5-^ - ((1 - VK)-^H,,9r,)) • det (l - VK)^,^^^ 

= (5-^ - ((1 - VK)-^Vfr„gr,)) ■ det (l - VK)^,^^^. 

(7.39) 

Since the second factor is independent of 5, the remaining task is the analytic 
continnation of the hrst. Using 115.5p . decompose f^ as 

1 r p-W^l3-riW^+sW 
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Now, 


(-^si 1) 9ri ) 




dZ 


f3+riZ^—sZ 

Z-6 


1 

27ri 



dZ 


^Z^/3+riZ^—siZ 

Z{Z-S) 


1 

6 


1 

27ri 



dZ 


^Z^l3+riZ'^-siZ 

Z{Z - 6) 


1 

6 


TZs- 


(7.41) 


The function TZs is analytic in 5 G M. Using these two identities as well as 
(1 — VK)~^ = 1 + (1 — VK)~^VK, we can rearrange the inner product as 
follows: 

= 1 _ ((1 _ VK)-^V1 + (1 - VK)-^Vf\gr,) 

= ^ - (PI + (1 - VK)-\VKV1 + Vr),9r,) 

= ^ - {Ps,l,9n) - {{V - Ps,)l + (1 - VK)-\VKV1 + Vn.9r,) 

= ns- ((1 - VK)-^ {vr + VKPs.l + (P - P.Jl), 

(7.42) 

Since is evidently analytic in 5 G M, we are left to show convergence of 
the scalar product. 

All involved functions are locally bounded, so to establish convergence it 
is enough to investigate their asymptotic behaviour, g^^ may grow exponen¬ 
tially at arbitrary high rate, depending on ri and 6, for both large positive 
and large negative arguments. We therefore need superexponential bounds 
on the function: 


(1 - VK)-^ {Pf* + VKPs,l + (P - P,J1). (7.43) 

For this purpose we first need an expansion of the operator P: 

n 

P = PiUri,r2 • • • -Psfe_iK-fc_i,rfcPsj;Urfc,ri- (7.44) 

k=l 

Notice that all operators P^., Pg. and Uri,rj map superexponentially decaying 
functions onto superexponentially decaying functions. Moreover Pg. and Pg. 
generate superexponential decay for large negative resp. positive arguments. 
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The function f* decays superexponentially for large arguments but may 
grow exponentially for small ones. Since every part of the sum contains one 
projection Pg^, Vf* decays superexponentially on both sides. 

Examining {V — Psi)!, notice that the k = 1 contribution in fl7.44p is 
equal to which is cancelled out here. All other contributions contain 
both Pg^ and P^^., which ensure superexponential decay. 

Using the usual asymptotic bound on the Airy function, we see that the 
operator K maps any function in its domain onto one which is decreasing 
superexponentially for large arguments. By previous arguments, functions 
in the image of VK decay on both sides, in particular VKPg^l. 

Now, in order to establish the hniteness of the scalar product, decompose 
the inverse operator as (1 —PK)~^ = 1 + VK[ \ —VK)~^. The contribution 
coming from the identity has just been settled. As inverse of a bounded 
operator, (1 — VK)~^ is also bounded. Because of the rapid decay, the 
functions P/*, VKPg^l and (P — Psjl are certainly in L^(M) and thus 
mapped onto L^(M) by this operator. Finally, the image of an L^(M)-function 
under the operator VK is decaying superexponentially on both sides. 

The expression 07.421) is thus an analytic function in 5 in the domain M. 
Setting (5 = 0 returns the value of Gm{T,V). Combining these results with 
07.39P hnishes the proposition. □ 

Lemma 7.5. The operator 1 — VK is invertible. 

Proof. We employ the same strategy as in |3]. For that purpose we use the 
following equivalence 

det(l + A) 7 ^ 0 1 + A is invertible. (7.45) 

Let 

( m 

Pl{A2(rfc) -rl< Sfc} 
k=l 

m \ 

Pl{A 2 (rfc) -rl< Smin} j > P(^max(A 2 (r) - r^) < 
k=l 7 

= pGOE(2^^^Smin) > 0 

(7.46) 

for any Smin > — cxo, where Pgoe is the GOE Tracy-Widom distribution 
function. For the last equality see [T2l[22]. The tails of the GOE Tracy- 
Widom distribution have been studied in great detail in various publications, 
see for instance 0. □ 
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8 Gaussian increments 


In this section we prove that the Airystat process has Brownian increments 
for nonnegative arguments: 


Theorem 8.1. Let 0 < ri < r 2 < • • • < Then 


P 


riM stat (rk) - A stat (rfc_i) e dak} 


k=2 


n 


^/ATT{rk - Tk-l) 


da. 


( 8 , 1 ) 


Proof. Without loss of generality we can assume that ri = 0. Denoting the 
partial derivative with respect to the Ath coordinate by di, we have 


P 


k=l 


riM stat ipk) ‘Pi j ^ ^ • • • 5 


i=l 


with 

A (si,..., Sm) = Gm{r, s) det (1 - VK) ^ 2 (k) • 
With a small abuse of notations, in what follows we will write 

( m \ / m 

n {-4 stat (rfc) e dsk] j = P riM stat {rk) = Sfc} 
k=l ' ^k=l 



Then, 


p riM stat (rk) = Sfc} )djA (si,...,, 
hfc=i / i=i j=i 

The crucial identity is: 


m m 


( 8 , 2 ) 


(8.3) 


(8.4) 


(8.5) 


IP f {"^stat (^fc) •^stat(^fe—l) ^k} j 

= f daiP^ n{Atat(’^fc) = CTI H-h ak}'] 

JK V / 

„ y m m \ 

= / dai ( TT ^ ) A (ai, ai + a2,..., ai H-h a„ 

Jm V / 

^ + - l-^m) 

m \ 

<9jjA(ai,ai + a2,...,ai + -- - + am) 


1=1 


cri=CO 


( 8 . 6 ) 


crj^ =—CO 
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We therefore have to study the asymptotics of A as Ui —)■ ±cxo. 
First we decompose A as 

A = Ai + As, 

Ai := (TZ — 1) det (1 — , 

As := det(l-VK- {Vf* + VKPs.l + (V - P.Jl) ® g) 

\ /L 2 (R) 


Since ri = 0 some functions simplify as 

poo poo 

n = Si+ dx dyAi{y), 

J Si j X 

/ oo 

dxki{x), 

poo ps 

g{s) = 1 — dxAi(a;) = / dxAi(a:), 

J S J —OO 

poo 

K{si,S2)= / dx Ai(si + a:)Ai(s3 + x), 

Jo 

where we used the identity (D.2) from |14j . 

Now consider Ai. 


( 8 . 8 ) 



(K-l)( det{l-VK) + din( det(l-PA:). 


2=1 


2=2 


(8.9) 

Regarding the hrst term, notice that the multiple derivative of the Fred¬ 
holm determinant gives exactly the multipoint density of the Airys process, 
which is known to decay exponentially for both large positive and negative 
arguments. This exponential decay dominates over the linear growth of TZ. 
Similarly, the {m — l)-fold derivative is smaller the (m — l)-point density of 
the Airys process, so this contribution vanishes in the limit, too. 

Continuing to As, using f* = —Kl, we hrst simplify the expression 


As = det — VK + {VKPg^l — (V — Psjl) <S) g^ ^ (8.10) 

We introduce the shift operator S, {Sf){x) = f{x + ai), which satishes 
SVr^,rjS~^ = and Pa+cri = S~^PaS, and consequently also 

1 “ Psi+(TiK-i,r2-Ps2+<Ti • • • Vrm-i,rmPsm+^i^r^,ri = ^ ^VS. (8.11) 


Using det(l — AB) = det(l — BA), we have 

As(s + cTi) = det (l - VSKS-^ + {VSKS-^Ps,l - {V - P,Jl) ®Sg^ ^ . 

( 8 . 12 ) 
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Now the dependence on the vector ^is only in the projection operators, while 
the dependence on cxi is only in these two operators: 


j-s+ux 

{Sg){s) = / da:Ai(x), 


' —OO 

poo 


(8.13) 


{SKS )(si,S 2 )= / dx Ai(si + a:)Ai(s 2 + a;). 


For large cxi, we have Sg ^ 1 and SKS ^ 0 (both strong types of con¬ 

vergence from the snperexponential Airy decay). So 

lim A 2 (s cTi) = det (l - (P - 0 l) = 1 - ((P - P^Jl, 1 )i, 2 (]k) 

0-1^00 V /L2(]g) 

(8.14) 

Applying the expansion 07.441) . we arrive at: 


n + CTi) = - ( JJ ) 5^(PsiKi,r2 • • • Ps^_xVr^_x,r^Ps^i, 1) 


2=1 


2=1 / k=2 


2=1 


( {Psi^rx,r2 ■ ■ ■ Psrn-l^r,n-l,rmPsm^i 1 ) 

(8.15) 

Writing ont this scalar prodnct and applying the fnndamental theorem of 
calculus leads to: 


(|8.15P V).i,r2(^l) ^2)h).2,r3 (^2) ■ ■ ■ ^rm-l,rmi^rn—ly ^m)y 


(8.16) 


which is the desired Gaussian density after setting Sj = X]I =2 as in 08.6p . 

For large negative cxi, we have Sg ^ 0 and SKS~^ —)■ 1. The rank one 
contribution is thus 

(PP,,1-(P-P,J1)0O. (8.17) 

We have to be somewhat careful here, as the convergence is weak (only 
pointwise) and (Sg) is not even L^-integrable. But the hrst factor decays 
super exponentially on both sides for hnite cxi and also in the limiting case 
PPs^l — (P — Psjl = (1 — P)Psil, so one should be able to derive nice 
convergence properties. Neglecting this rank one contribution we are left 
with 

= 0. (8.18) 

□ 


lim A 2 (s-|- (Ti) = det (l — Pi) 

!->■ —OO V / 


L2(]R) 
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